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Abstract. We propose a generalized framework for the study of voter models in complex networks at the 
heterogeneous mean-field (HMF) level that (i) yields a unified picture for existing copy/invasion processes 
and (ii) allows for the introduction of further heterogeneity through degree-selectivity rules. In the con- 
text of the HMF approximation, our model is capable of providing straightforward estimates for central 
quantities such as the exit probability and the consensus/fixation time, based on the statistical properties 
of the complex network alone. The HMF approach has the advantage of being readily applicable also in 
those cases in which exact solutions are difficult to work out. Finally, the unified formalism allows one to 
understand previously proposed voter-like processes as simple limits of the generalized model. 
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1 Introduction 

A topical problem in the statistical physics approach to 
social and evolutionary dynamics |1I2| is the study of the 
mechanisms ruling the formation of consensus in an ini- 
tially disordered population, in situations implying the 
opinion about a certain issue, the intention of voting in 
an election, or the evolutionary competition of different 
species striving for the same ecological resources. Several 
stochastic copying/invasion processes have been proposed 
to represent this kind of problems, the simplest being the 
voter model [3] and the Moran process [4]. In these models, 
each individual in a population (agent) is endowed with 
a binary variable (opinion or state) with value a = ±1. 
At each time step, an agent i, together with one nearest 
neighbor j, are selected at random. In the voter model 
the system is updated as Oi :— cr,-, the first agent copying 
the opinion of its neighbor. The Moran process, on the 
other hand, can be considered as a reversed voter model, 
in which it is the neighboring agent the one who copies 
the opinion of the first agent, Oj :— <jj. Starting from 
a disordered initial state, this kind of dynamics leads in 
finite systems to a uniform state with all individuals shar- 
ing the same state (the so-called consensus). The way in 
which this final state is reached is usually characterized 



in terms of the exit probability E(x) and the consensus 
time Tn(x), defined as the probability that the final state 
corresponds to all agents in the state +1 and the aver- 
age time needed to reach consensus in a system of size 
N, respectively, when starting from a homogeneous initial 
condition with a fraction x of agents in state +1 [T]. 

While simple voter-like models are well understood 
on regular lattices, even in terms of theorems and exact 
solutions |5l6j . they become more relevant in social and 
evolutionary contexts when considered on top of complex 
networks, which act as more realistic representations of 
social or ecological contact patterns [7|8I9] . The analy- 
sis of the voter model in these substrates reveals nontriv- 
ial differences with respect to ordered lattices. For exam- 
ple, now the order in which interacting individuals are 
selected becomes relevant [10], in such a way that voter 
model and Moran process behave in different ways. More- 
over, relevant quantities such as the consensus time turn 
out to depend on the heterogeneity of the contact pat- 
tern, as measured by the degree distribution |11I12I13] , 
Not only the properties of the interaction substrate are 
relevant in this case, but also the intrinsic heterogene- 
ity of the actors. Their individual propensity to interact 
with peers, and change state accordingly, plays a signifi- 
cant role [14115116117] . Other sources of heterogeneity that 
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have been considered include explicitly the directed nature 
of the connections [18] . 

The theoretical understanding of voter-like models (and 
dynamical processes in general) on complex networks has 
been traditionally accomplished by application of hetero- 
geneous mean- field (HMF) approaches [19120) . which are 
based on a twofold assumption: (?) The network descrip- 
tion is coarse-grained into degree classes, all vertices in the 
same class having the same degree and sharing the same 
dynamical properties; (ii) The real (quenched) network 
structure is replaced by an annealed one |20) , which disre- 
gards the actual connection pattern and simply assumes 
that the degree class k is connected to the degree class 
k' with conditional probability P(k'\k) [21 . Very signifi- 
cant progress has been achieved in the analysis of voter- 
like models within the HMF approach, which allows one 
to work out simple analytic expressions for the quantities 
of interest, showing reasonable agreement with numerical 
simulations in real quenched networks [11112113114117122"] . 

Recently, a generalized formalism for the class of het- 
erogeneous stochastic-copying voter-like models on net- 
works has been proposed [23], in which the process is 
identified by the copying rate Cy , encoding the full struc- 
ture of the contact network and the stochastic update 
rules, and defined as the rate at which vertex i in the 
network copies the state of vertex j. Thus, for exam- 
ple, the standard voter model corresponds to the choice 
Cij — aij/[Nki] with ki = 2^ r a ir; where ay is the adja- 
cency matrix of the network and N the network size. The 
Moran process, analogously, is given by Cy = aij/[Nkj]. 
Within this formalism, it has been shown that both the 
exit probability and the consensus time can be calcu- 
lated exactly from the knowledge of the spectral prop- 
erties of the matrix Cy, provided that certain general 
conditions are met [23124] . a result that lays the foun- 
dation for a mathematical understanding of general copy- 
ing processes and their mapping to particle-reaction sys- 
tems. Despite the fact that the formalism in Refs. [23125] 
is exact, and provides in some cases more accurate re- 
sults than HMF theory, it is still useful to consider general 
stochastic-copying models from the perspective of HMF 
theory. In this framework, indeed, approximate analyti- 
cal results can be obtained when the exact solution would 
be hard to work out in practice. For example, in realis- 
tic heterogeneous environments involving large numbers 
of agents, explicit expressions for Cy might not be readily 
accessible. Moreover, the spectral properties of the copy- 
ing rates are in general non-trivial to obtain, unless the 
matrix Cy has a relatively simple form. 

In this paper we pursue this path, proposing a gen- 
eralized coarse-grained voter-like model on networks and 
showing how the HMF approach allows us to obtain very 
simple estimates for central properties such as the exit 
probability and the consensus time. We check the valid- 
ity of our approach by considering a simple example of 
opinion dynamics in a homophilic society, in which ver- 
tices with similar degree are more prone to interact than 
vertices with differing degree. 



2 Generalized voter model on networks 

Inspired by Ref. [23], we consider a stochastic model on 
networks defined in terms of a heterogeneous voter model 
as follows: 

— Each vertex i is endowed with a given fitness /, [12] . 

— A source vertex i is selected at random, with a proba- 
bility fi/J2j fji i- e -j proportional to its fitness fi. 

— A nearest neighbor j of i is then selected at random. 

— With probability Qy, i copies the state of vertex j. 
Otherwise, nothing happens. 

With these settings, the microscopic copying rate Cy, as 
considered in Ref. [23 will be given by 



Cii — 



fi 



J2j fj k-< 



-Qi 



(i) 



In the spirit of the HMF approximation, we can replace 
the microscopic copying rate by its degree class average. 
The quantities fi and Qy are simply coarse-grained by 
averaging them over the set of vertices with a given fixed 
degree, i.e. 



fi 

Qij 



1 



NP{k) 



} ] fi — fk i 



iefc 



1 



NP(k) NP(k 



(2) 

knEE^-^^o, ( 3 ) 



where i £ k denotes a sum over the degree class k and 
P(k) is the network's degree distribution. For the term 
concerning the random choice of a nearest neighbor, we 
follow Ref. 126 to substitute 



ki 



— > 



[NP(k)]-'j: iek z 



jek' ai i 



[NP(k)]-^iekEr^r 



P{k'\k). (4) 



At the coarse-grained degree level, our generalized voter 
model is thus defined in terms of the mesoscopic copying 
rate 



C(k,k') 



/(fc) 
</(fc)> 



P(k'\k)Q(k,k'), 



(5) 



where the function f(k) comes with its proper normaliza- 
tion factor. 



3 Heterogeneous mean-field solution 

In the HMF approach, ordering processes are quantified by 
studying the evolution of the density of vertices of degree k 
in the state +1, Xk- In order to determine the rate equation 
satisfied by these quantities [11113117] . we consider the 
probability i7(fc; a) that a spin in state a at a vertex of 
degree k flips its value to — er in a microscopic time step. 
From the definition of the generalized voter model, this 
probabilities can be simply written as 



n(k;a) 



I- a 



E 



i 



ax k > 



P(k)C(k,k'). 
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From the previous expression, the rate equation for x k can 
be written as IT7l 



Xk 



n(k;-i)-n(k-,+i) 

P(k) 



The model just posed is still too hard to solve even in 
the HMF approximation. In the following, we will there- 
fore make two major simplifying assumptions leading to 
analytically solvable HMF equations: (i) Dynamics pro- 
ceed on uncorrelated networks, i.e. [5] 



As for the usual voter model the conservation law 
allows the immediate determination of the exit probability 
E, i.e. the probability that the final state corresponds to 
all spins in the state +1. In the final state with all +1 
spins we have uj = 1, while ui = is the other possible 
final state (all —1 spins). Conservation of ui implies then 
w = E ■ 1 + [1 - E] ■ 0, hence 



P(k'\k) 



k'P(k') 

(k) ! 



(7) 



and (ii) the interaction probability can be factorized as 

Q{k,k l ) = a{k)b{k')s{k : k'), (8) 

where s(k,k') is any symmetric function of k and k' . As 
we will see, this simplified form, which permits an analytic 
HMF solution, still allows for a vastly rich phenomenology, 
encompassing all voter-like models on complex networks 
previously proposed. Under this conditions, defining 



u{k) 



a(k)f(k) 



v(k') 



b(k')k' 



(/(*)> ' " V " J (k) 
the HMF rate equation can be written as 



■J'k 



= J2P(k')r(k,k')( 



Xk) 



(9) 



(10) 



where r(k, k') = u(k)v(k')s(k, k'), 

The HMF analysis proceeds by first determining the 
corresponding conservation laws [11112113] . Conserved quan- 
tities for the generalized process in Eq. (fTU)) can be cal- 
culated as follows: We define a generic integral of motion 
w[xfe(t)] such that ckv/dt = 0. By definition of time deriva- 
tive, we have 



duj 



Edtu . 

k 



0. 
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In analogy with previous results [11112113] . we look for 
conserved quantities that are linear in Xk imposing dui/dxk 
Zk independent of x k , so that conserved quantities will be 
given by 

" (12) 



E = uj = 



(v{k)/u(k) x k ) 
(v(k)/u(k)) 



(14) 



Starting from a homogeneous initial condition, with a given 
density x of randomly chosen vertices in the state +1, we 
obtain, since ui = x, 



E h {x) 



(15) 



completely independent of the defining functions a, b, and 
s, and taking the same form as the standard voter model 
[I]. On the other hand, with initial conditions consisting 
of a single +1 spin in a vertex of degree k, we have 



E x {k) = 



v(k)/u(k) 
N(v(k)/u(k))' 



(16) 



which does not depend on the functional form of the sym- 
metric interaction term s(k, k'). 

By looking at Eq. ([TO]) , every choice of Xk constant in k 
is a solution to the steady state condition x k = 0. We can 
prove that this solution is unique and does not depend 
on initial conditions if the square matrix P{k')r{k, k') 
is irreducible and primitive (it certainly is when work- 
ing with positive rates, which we will do in the follow- 
ing) [27]. We shall call the solution for the steady state 
Xk{t — >• oo) = x 00 . Then it is easy to prove that 



E 



Zk'Xk' 



(17) 



Z • X 



E 

k 



ZkXk ) 



where Zk is any solution of ^2 k Zkik — and ik is given 
by Eq. (JTUJ). Considering the explicit form of Eq. (fTU]) . 
the choice Zk oc P(k)v(k)/u(k) always satisfies the above 
condition, so that a conserved quantity is found up to 
multiplicative factors and additive constants. We choose 
the normalization z k = 1, such that the conserved 
quantity is defined as 



that is, even in this general case, the steady state value 
for Xk equals the conserved quantity. 

After a rapid exponential convergence to the steady 
state distribution, the systems starts fluctuating diffusively 
around this value, until consensus is reached. Such fluctu- 
ations characterize finite systems and occur at long time 
scales, making such two-step relaxation process possible 
in most cases. The average consensus time Tjv(x) for a 
system in a generic state x can be derived extending the 
well known recursive method to our general case [13] . At 
a given time t, Tjv(x) must equal the average consensus 
time at time t + At plus the elapsed time At = 1 /N that 
is, in our notation, 



Tjv(x) =nT N (x) 



-Y,n(k;s)T N ( X - 

k,s 



Ax (fe) )+Z\i, (18) 



(v{k)/u{k)x k ) 
(v(k)/u(k)) 



(13) 



where 77 = 1 — J2k s s) is the probability that no 
state change occurs, while the sum is the weighted average 
over possible state-updates x — > x+ Ax' 4 '. The variation 
Ax'*' is a vector whose all components are zero except for 
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Model 


/(*) 


a(k) 


b(k) 


s(k,k') 


Voter model [3] 


1 


1 


1 


1 


Moran process [4] 


1 


k 




1 


Link update [10] 


k 


1 


1 


1 


Voter weighted [T7] 


1 


1 


9a(fc)(fc> 
(fe 9s (fc)> 


1 


Moran weighted [T7] 


1 


kg s (k) 
(fcg s (fc)> 


(fc) 

k 


1 


Generalized Voter [TJ] 


1 


1 


k"- 1 


k + k' 
k a +k'<* 



Table 1. Summary of the mapping of voter-like models to the 
present formalism. 



the fc-th, which equals the update-unit A k = [NP(k)}^ 1 . 
Expanding to second order in A k , taking x k = lu as the 
initial state and changing variables such that d/dxk = 
z k d/duj we obtain the backward Kolmogorov equation 



1 = 



z J Tz 



u(l — uj) 



d 2 r N 



N ~ v ~ "' dtu 2 
leading to 

Tn(uj) = — N c h[uj hi lj + (1 — u>) ln(l — u>) 



(19) 



(20) 



where we have defined the effective system size iVeff = 
N/J2kk' z k.r(k,k')zk> , which, in the case of generalized 
voter dynamics, Eq. ([9]), becomes 



N 



(f(k))(k) 



kb(k) \ A 
/(fc)o(fc)/ 



s(k, k')kb(k) j (fc /) a( fc, 



(21) 



with ((•))= J2 kk ,P(k)P(k%). 

Equations (|2T)]) and (|2"Tj) . together with the expres- 
sion for the exit probability, Eq. (jl~4"]) . represent the fi- 
nal HMF solution of the generalized voter model. From 
these formulas it is easy to recover most of the varia- 
tions of the voter model considered in the past. For exam- 
ple, the standard voter model is obviously recovered for 
a(k) = b(k) = /(fc) = ,s(fc,fc') = 1. The invasion process 
[28 12 , also known as the Moran process in the evolution- 
ary literature |4l29j . corresponds to a(k) = fc, b(k) = fc -1 
and /(fc) = s(k, fc') = 1. Link update dynamics [10] is re- 
covered for a(k) = b(k) = s(fc, fc') = 1 and /(fc) = fc. The 
voter and Moran processes on weighted networks charac- 
terized by a symmetric weight between vertices of degree fc 
and fc' proportional to g s {k)g s (k') [T7] are reproduced by 
imposing s(k, fc') = 1 and setting a(fc) = /(fc) = 1, b(k) = 
g s (k)(k)/(kg s (k)) and /(fc) = 1, a(fc) = kg s (k)/{kg s (k)), 
b(k) = (fc)/fc, respectively. Finally, an HMF implementa- 
tion of the generalized voter dynamics proposed in [14] is 
recovered imposing /(fc) = a(k) — 1, b(k) — fc Q_1 and 
s(fc, fc') = (fc + fc')/(fc Q + k' a ). A summary of the above 
mappings is presented in Table [1] 



4 Numerical analysis 

In order to show an application of our formalism, we ex- 
amine a toy model to study the effects that homophily 



in social networks might have in opinion formation dy- 
namics. Homophily is broadly defined by the tendency of 
people to interact with similar people |30| . In the absence 
of information beyond the topological structure of the con- 
tact network, the simplest assumption we can make is that 
homophily is driven by an increased tendency of individ- 
uals to copy other individuals with a degree that is not 
too different from their own. As a simple representation 
of homophilic behavior we consider the case s(k, fc') = 
exp[-i?(fc - fc')/£ 2 ] and /(fc) = a(fc) = b{k) = 1. Here 
R(x) is any continuous even function with R(0) = and 
a minimum in x = 0, i.e. R'(Q) = and R"(Q) > 0. The 
parameter £ measures the amplitude of stochastic fluctua- 
tions around ideal homophily: For £ w 1 the probability of 
copying neighbors with different degrees is strongly sup- 
pressed; for £ 3> 1 fluctuations take over and simple voter 
behavior is rapidly recovered. 

In spite of its apparent simplicity, this problem would 
be impossible to solve in a realistic network by standard 
techniques [53]. By applying the HMF result in Eq. ([2"T]) 
instead, we readily find 



7V« = N- 

eft 



EfeEfe' fcP(fc)fc' 2 P(fc')exp 



R(k-k') 



(22) 



where we consider complex networks with a scale-free de- 
gree distribution P(k) ~ fc~ 7 , with 7 s]2,3], minimum 
degree m and maximum degree fc c . In the limit of small 
£ the denominator of Eq. (|2"2")l can be evaluated in the 
continuous-degree approximation. By applying the Laplace 
method for the first term in the asymptotic expansion, af- 
ter an adequate change of variables, it is thus straightfor- 
ward to find that for small £ 



Ni a ~ — 



N 
I 
2(1 



R"(0) 



2tt 

7) 



(k 



2-7 



(2 - 7 ) 2 [kl^ - m i-r)(k 2 c -~> + m 2 -y) 



(23) 



whereas in the opposite limit £ — >■ 00 one recovers the 
simple voter model result, which in our notations reads 
N^g ~ Nk2~ 3 - As a consequence, the consensus time 

Tn oc N^ s diverges for small £ as £ -1 as the selectivity 
amplitude £ approaches zero. For increasing £, instead, 
Tn decreases and asymptotically crosses over to a plateau, 
where simple voter behavior is recovered. 

In order to check explicitly the predictions of our for- 
malism, we consider a Gaussian homophily model, given 
by the simplest choice R(x) = x 2 . In Fig. [5] we plot the 
consensus times for homogeneous initial conditions as a 
function of £, computed from the numerical evaluation of 
Eq. (|2"2"]) . We observe that, starting from large values of 
£, Tat is constant, as expected for simple voter behavior. 
Upon decreasing the selectivity amplitude £, the consensus 
time starts increasing, asymptotically behaving as As 
soon as £ decreases below one, the discrete nature of the 
degree distribution takes over and TV reaches a plateau. 
Interestingly, we find that TV is an increasing function of 
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for large £. At intermediate values of £, a reasonable fit 
to the form can be obtained. For smaller values of 
£, the selectivity range is so narrow that quenched effects 
of the network topology take over, tending to slow down 
the dynamics beyond the HMF prediction [T7]. Values of 
Tjv fluctuate wildly and deviate from the predictions for 
annealed-network topologies. 



5 Conclusion 



10 



10 



10 
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Fig. 1. Consensus time as a function of the selectivity range 
from numerical evaluation of Eq. (|22|l . The case of homoge- 
neous initial conditions Xfc(0) = 1/2 (maximum entropy) is 
considered. Inset: Consensus time as a function of the degree 
distribution exponent. Numerical estimates (points) and ana- 
lytical predictions (lines) based on N£g and AT 
text are compared. 
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Fig. 2. Consensus time with a Gaussian homophily factor in 
scale-free networks with 7 = 2.5, m = 4 and k c ~ N 1 ^ 2 . Main 
plot: Consensus time as a function of N, for different values of 
£. Inset: Consensus time at constant N = 5000 as a function 
of £. The full line is a guide to the eye with slope —1. 



7 for large £, while it decreases with 7 in the small £ limit 
(Fig.Hl inset). 

These observations have also been checked against di- 
rect numerical simulations of the Gaussian homophily mo- 
del on uncorrelated scale-free networks, generated with the 
uncorrelated configuration model [21] • In Fig. [5] we show 
(main plot) the consensus time as a function of the net- 
work size N for different values of £. We can see that the 
overall plot increases at fixed N for decreasing values of £, 
while for large £ it tends to the limit given by the simple 
voter model. When considering Tjy as a function of £ in 
a network of fixed size, inset of Fig. we observe, as ex- 
pected, an increase for decreasing £, followed by a plateau 



In this paper we have studied a generalized voter-like 
model in the framework of the HMF theory. Remarkably, 
the HMF formalism we have adopted allows for a straight- 
forward mapping of most of the voter-like models proposed 
in the past. Not only the the properly said voter model, 
the Moran process, and their weighted generalization can 
be reproduced in our formalism, but also the link-update 
dynamics and recently introduced generalizations of the 
voter model can be straightforwardly recovered. The HMF 
approach allows for predictions that are in fair agreement 
with the exact results whenever they are available, and 
most importantly provides approximate solutions in the 
event that exact methods are not viable. 

As an instance of the versatility of the HMF approach, 
we have considered a simple example of a heterogeneous 
voter model in which vertices with similar degree tend to 
interact more often among themselves rather than with 
the rest of the network. This degree selectivity is a natu- 
ral HMF implementation of the concept of homophily, i.e. 
the tendency of individuals to associate and bond with 
similar others. While an exact solution to such problem 
would be unworkable, we have shown that in this frame- 
work homophily ends to slow down consensus dynamics, 
in reasonable agreement with the predictions of HMF the- 
ory. 

In the future, it would be interesting to extend the 
approach to voter-like methods that exhibit a surface ten- 
sion, as for example the Naming Game 32 , or the noise- 
reduced voter model |33| . as well as to improve the qual- 
ity of our present results by applying more sophisticated 
HMF approaches, such as the dynamical pair approxima- 
tion |34) . or the master equation approach |35) . 
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